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We derive and analyze the Born-Markov master equation for a quantum harmonic oscillator 
interacting with a bath of independent two-level systems. This hitherto virtually unexplored model 
plays a fundamental role as one of the four "canonical" system-environment models for decoherence 
and dissipation. To investigate the influence of further couplings of the environmental spins to 
a dissipative bath, we also derive the master equation for a harmonic oscillator interacting with 
a single spin coupled to a bosonic bath. Our models are experimentally motivated by quantum- 
electromechanical systems and micron-scale ion traps. Decoherence and dissipation rates are found 
to exhibit temperature dependencies significantly different from those in quantum Brownian motion. 
In particular, the systematic dissipation rate for the central oscillator decreases with increasing 
temperature and goes to zero at zero temperature, but there also exists a temperature-independent 
momentum-diffusion (heating) rate. 



PACS numbers: 03.65.Yz, 42.50.Lc, 85.85.+J, 37.10.Ty 



I. INTRODUCTION 

Theoretical studies of decoherence and dissipation in 
quantum systems have hitherto focused on three canoni- 
cal system-environment models: (i) A harmonic oscilla- 
tor (or, more generally, a particle moving in phase space) 
coupled to a bath of other harmonic oscillators (quantum 
Brownian motion) [l|, 0, y, |j| ; (ii) a quantum two-level 
system (TLS), represented by a spin-i particle, inter- 
acting with a bath of harmonic oscillators (spin-boson 
model) Q; and (iii) a spin-i particle coupled to a bath 
of other spins (spin-spin model) [fj. Surprisingly, how- 
ever, the fourth possible canonical combination, namely, 
a single harmonic oscillator interacting with a bath of 
spin-i particles — which, in obvious nomenclature, shall 
henceforth be referred to as the oscillator-spin model 
(Fig. QJi) — has not yet been studied in any detail in the 
literature in terms of a Markovian master equation. 

It is the purpose of this paper to close this gap by giv- 
ing a microscopic treatment of the oscillator-spin model. 
We will derive the Born-Markov master equation and 
compare the resulting dynamics to those induced by an 
oscillator bath (quantum Brownian motion) . Apart from 
its relevance in completing the set of canonical mod- 
els, the oscillator-spin model is also motivated by re- 
cent experiments on quantum-electromechanical systems 
(QEMS) 0, S i and micron-scale ion traps. In both 
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systems, a central quantum-mechanical vibrational de- 
gree of freedom interacts with two-level defects causing 
dissipation and decoherence of the oscillator. We may 
represent this situation by a harmonic oscillator coupled 
to a collection of TLS, i.e., by a model of the oscillator- 
spin type. 

QEMS are nanometer-to-micrometer-sized crystalline 
mechanical resonators coupled to nanoscale electronic 
transducers that detect the high-frequency (MHz-GHz) 
vibrational motion of the resonator. Since only the low- 
est, fundamental flexural mode of the resonator turns 
out to be relevant Q, the resonator can be modeled as a 
single quantum-mechanical harmonic oscillator. Recent 
experimental evidence 1C , [ll[ (see also the molecular- 
dynamics simulation of 12 1 and earlier results in p^LH3|) 
strongly suggests that the dominant source of decoher- 
ence and dissipation in QEMS is the interaction with 
two-level defects intrinsically present in the resonator it- 
self. 

In ion traps a single ion can be trapped by a time- 
dependent potential and cooled to very low energies [l5[ . 
Under the right conditions the motion of the ion is equiv- 
alent to that of a quasi-one-dimensional harmonically 
bound particle. A major source of decoherence in ion 
traps is thought to arise from fluctuating patch poten- 
tials on the trap electrodes [l6[. Roughly speaking, this 
causes a fluctuating linear potential that results in ran- 
dom forces acting on the ion. The net effect is a slow 
heating of the ions. This problem is particularly acute in 
small traps [13, EU , where anomalous heating has been 
experimentally observed fl9j | . More recent experiments 
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of models tailored to specific experimental situations. In 
existing models of ion traps, the fluctuatin g fo rces have 
thus far been treated classically IM III III SI ■ For 
QEMS, a realistic and quantitatively accurate model- 
ing of the influence of the various defects on the res- 
onator is rather involved. First theoretical studies (see, 
e -g-; [Hi EH) were recently followed by detailed work 
by Seoanez, Guinea, and Castro Neto [27|, [28|. In order 
to be able to use a spectral-function treatment for the 
environmental TLS, these authors focused on the limit 
fceT > where tto is the natural frequency of the 

resonator. However, attaining the quantum regime of 
QEMS requires the opposite limit k^T <C hfl (GHz 
QEMS). In this case, it is likely that only very few TLS 
will be relevant to the dynamics of the resonator and we 
suggest that the correct description may be closer to the 
model discussed in Sec. Mil below. 

This paper is organized as follows. Sec. |TT] presents 
the derivation of the Born-Markov master equation for 
the central oscillator coupled to a spin bath. In Sec. ITTT1 
we derive the master equation for a harmonic oscillator 
interacting with a single spin coupled to a bosonic bath. 
We summarize our results in Sec. IIVI 



FIG. 1: (a) Oscillator-spin model, (b) Oscillator-spin model 
with each spin coupled to an additional bosonic bath, (c) 
Limiting case of (b) in which the oscillator interacts with a 
single environmental spin coupled to a bosonic bath. 



II. MASTER EQUATION FOR A HARMONIC 
OSCILLATOR COUPLED TO A SPIN BATH 

A. Model 



have cooled the traps to a few K and seen a dramatic 
reduction in heating [20(. We anticipate that with fur- 
ther cooling the heating will ultimate be attributable to 
charge fluctuations in two-level traps, especially for oxide 
barriers in semiconductor substrates. 

In many cases, it is reasonable to assume that each 
of the two-level defects will be also coupled to its own 
environment, which we may model as a bosonic bath. 
We are thus led to a more complicated model in which 
the central oscillator couples to a collection of indepen- 
dent spin-boson models (Fig. [Tp). In this paper we will 
consider the special case of only a single TLS interact- 
ing with a bosonic bath (Fig. [T}:). This simplification 
allows us to analytically derive the master equation for 
the central oscillator in the limit in which the oscillator is 
strongly coupled to a steady-state TLS. It is also exper- 
imentally motivated by the fact that in GHz QEMS and 
micron-scale, cryogenic ion traps the number of defects 
that participate in the dynamics is thought to be quite 
small. 

We emphasize that the focus of this paper is a study 
of the general dynamics and properties of the oscillator- 
spin model in the context of the canonical models, and 
it is not our aim to present detailed models for ion traps 
and QEMS. However, these systems lend urgent exper- 
imental relevance to the oscillator-spin model, and our 
model may serve as a starting point for the development 



We consider a single quantum harmonic oscillator (the 
system S) with sclf-Hamiltonian 
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(1) 



The oscillator interacts with an environment £ of N in- 
dependent spin-^ particles (quantum TLS). The environ- 
ment is described by the self-Hamiltonian (setting h = 1) 

N N N . 

^^^-E^'+Et^ ( 2 ) 



i=i 



i=l 



where u>i and A, are, respectively, the asymmetry energy 
and tunneling matrix element of the ith bath spin. The 
environment couples linearly to the position coordinate 
of the oscillator via the interaction Hamiltonian 



H int = X®J2^i l) = X®E. 



(3) 



i=l 



The total system-environment combination is then de- 
scribed by the Hamiltonian 



H — Hs + Hf + H\, 



(4) 



We assume the limit of weak system-environment cou- 
plings and no initial system-environment correlations, 
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p(Q) = ps(0) <8> /5g(0). We take the environment to be in 
thermal equilibrium at temperature T. Since the spins 
of this thermal bath are independent, we have 
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i=l 



(5) 



where j3 = 1/fceT and Z = Ti£ e^ l3He . We would now 
like to derive the Born-Markov master equation for this 
spin-bath model. For the interaction Hamiltonian ([3]), 
the general form of the master equation is [29j 



d7^ 



H s ,ps(t) 

&t!c(t) X,X(- T )p s (t) 



+ C(-t) p s (t)X(-r),X 



(6) 



Here 



X(t)=X cos (fi r) + — - P sin (fi r) (7) 

denotes the system's position operator X in the interac- 
tion picture. The spin-environment self-correlation func- 
tion C(t) appearing in Eq. (J6j) is given by 



C(r) = 



where -E(r) = e lH£T Ee lH£T is the environment operator 
_E in the interaction picture and the average is taken over 
the initial state pg = ps{0) of the environment (the Born 
approximation means that p(t) « ps(t) <8 pe (0) for all t). 



B. Calculation of the environment self-correlation 
function 

First, we compute the environment self-correlation 
function |5]). which we may write as 



Pe 



U) 



(9) 



Because the environmental spins do not directly interact 
with each other, they are uncorrelated, 



(<7«(r)a«)_ =(<t«(t)\ for j, 

\ I p £ \ ' pe ^ ' pe 



and thus Eq. © can be decomposed as 



(10) 



P£ . x 1 Pe 

3^ 



Let us assume that at t = the average of the "quantum 
force" due to the collective action of all environmental 
spins is equal to zero, 



(12) 



' pe x ' pe 



This is a nonrestrictive assumption, since Eq. (112|) can 
always be fulfilled by simply adding a constant to the 

Hamiltonian. Then the term 5j \&z ) appearing 

\ I pe 

in Eq. (fTTj) will also tend to zero, and Eq. ([9]) simplifies 
to 



C(r)=Es^f'(^ 



i JJe-^ <)/fcBT 



« 



= £rfiM e ~** )/teT ^°< T ^ fl }' (13) 

i 

where Z$ = Tr^ e~ H£ «' feBT anc j i n the last line we have 
again used the fact that the bath spins are uncorrelated. 

To calculate the interaction-picture operator ai (t) = 
e -iH e fjW e -'^ £ ^ we wr jte the environment Hamilto- 
(8) nian Hg in matrix form in the eigenbasis {\0) i: |l)j} of 
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(14) 



The matrix eigenvalues are E±' = ±^^/u>f + A| = ±5^ 
with corresponding eigenvectors 

l+) ! -cos||0) l + sin||l) 4 , (15a) 
R^-sinllO^+cosIll),, (15b) 

where 0t = arctan ^ . With respect to the basis 
I - )i) ) the matrix representation of eri 1 ' reads 



. (i) / cos Pi — sin ( 
- sin f i — cos ( 



(16) 



We can now evaluate Eq. (113|) directly by carrying out 
the relevant matrix products and then taking the trace. 
The result is 



C(r)=C + J2 



[COS (uJiT) 



-itanh(/35 l /2)sin(w l r)], (17) 
where C = J2i i^ff^ i s a time-independent constant. 
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Continuum limit 



Here, the coefficients f^, 7, D, and / are defined as 



Let us introduce the spectral density function 



m = E 

i 

and write Eq. IT7|) as 



5(lo - u>i) 



(18) 



C(r) = C + I dui J (to) [cos (lot) 





— i tanh(/3u> /2) sin (wr)] 
= Co + i/(t) - hj(t). (19) 

The functions i/(r) and 7y(r) take the same functional 
form as the noise and dissipation kernels, respectively, in 
the case of an oscillator bath (quantum Brownian mo- 
tion) with spectral density J osc (lo), 

^osc — J dco J osc (lo) coth t ) C ° S ^ T ' ) ' '' 20a '' 



Vosc(t) = I dco J osc (ui) sin (lot) , 
10 



provided we choose 



J (no) = Jqsc(w) coth 



2k B T 



(20b) 



(21) 



in Eq. (|19p . Conversely, ignoring the constant term Co, 
we can map the spin bath with spectral density J(lo) onto 
an oscillator bath with "surrogate" spectral density 



Josc(lo) = J '(to) tanh 



2k B T 



(22) 



This is an example of the general result, first derived 
by Feynman and Vernon [30(, that in the limit of suf- 
ficiently weak coupling any dissipative bath (including 
the spin bath) can be mapped onto a bath of oscillators. 
We note that in the limit Aj 3> u>i for all i (and thus 
Co — » 0), expression (fT9)) coincides with a result previ- 
ously obtained by Caldeira, Castro Neto, and de Car- 
valho [31] for a model of a general system interacting 
linearly and weakly with a spin bath, where the envi- 
ronmental self-Hamiltonian was assumed to take a more 
simple form than in our model. 



D. Master equation 

Inserting Eq. (fn?|) into Eq. © and using Eq. (O leads 
to the master equation 



d 
dl 



Ps(t) 



H s + hm 2 Q x 2 ,p s (t) 
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,{P,Ps(t)} 
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X, 



X,p s (t) 



f X, P,p s (t) 



(23) 



fi 2 
"0 



"M J dr? ?( T ) cos (^oT) 



dr tj(t) sin (ilo r ) 



(24a) 
(24b) 



D= dr [C + v(t)] cos (n r) 
Jo 

poo 

= Dq+ dr v(t) cos (Qqt) 



D, 







JO 
Di, 



(24c) 



1 



MQ j 
f 1 

/0 + /1. 



dr [C + f(r)] sin (fi r) 
dr !/(r) sin (Qqt) 



(24d) 

The interpretation of these coefficients is analogous to 
the case of quantum Brownian motion. The coefficient 
Slg describes a frequency shift ("Lamb-shift" renormal- 
ization of the natural frequency of the oscillator), 7 is the 
momentum-damping (and thus dissipation) rate, and D 
and / are the normal-diffusion and anomalous-diffusion 
coefficients describing decoherence. We see that Sl§ and 
7 are explicitly temperature-dependent while D and / 
are not, which is exactly opposite as in quantum Brow- 
nian motion. Formally, this difference is easily under- 
stood from the fact that using the surrogate spectral den- 
sity (12"2")) in the expressions for the oscillator-bath noise 
and dissipation kernels (|20p eliminates the temperature- 
dependent term in the integral (|20ap while introducing 
the term tanh (fl /2k B T) in the integral (|20b[) . 



E. Example: Ohmic spectral density 

Let us consider an ohmic spectral density for the spin 
bath with a Lorentz-Drude high-frequency cutoff, 

A 2 



2Af 7o ^ 

J(U0) = LO 



A 2 



-a ■ 



(25) 



The coefficient 7 is given by a double Fourier sine trans- 
form of the function J (no) tanh(/3a5/2), which returns the 
original function up to a prefactor of tt/2, 



7 



1 f°° 

— dr sin (O r) 
"0 Jo 



x / duj J(lo) tanh(/3w/2) sin (ujt) 
Jo 

Similarly, the coefficient D\ is given by a double Fourier 
cosine transform of J (to) , which leads to 



D = D + Mj n 



A 2 



A 2 + ft 2 ' 



(27) 
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FIG. 2: Dissipation and normal-diffusion coefficients (a) 7 
and (b) D as a function of the bath temperature T for the 
oscillator-spin model (solid line) and for quantum Brown- 
ian motion (circles), assuming the ohmic spectral density 
(|25[) . The vertical axis is normalized in units of the zero- 
temperature values j(T = 0) and -Dqbm (T = 0), respec- 
tively. The horizontal (temperature) axis is displayed in units 
of n /k B - We use D = Dqbm(T = 0). 



For quantum Brownian motion these coefficients read [2£ 



A 2 



7QBM = 70 



A 2 + fi§' 



-Dqbm = Mj tl 



— ■= ^ coth | — ) 

A 2 + fl 2 \2k B TJ 



(28a) 
(28b) 



Disregarding the constant term Do, we see that the spin- 
bath coefficients are given by the oscillator-bath coeffi- 
cients multiplied by the term tanh (Qo/2k-BT), which is 
simply a direct consequence of the use of the surrogate 
spectral density 



temperature. This initially surprising result is easily ex- 
plained. While a harmonic oscillator can absorb an infi- 
nite amount of energy, there are only two energy levels for 
a spin-l particle. It follows that, as the bath temperature 
is raised, the spin bath saturates quickly and the dissipa- 
tive influence on the central system must decrease when 
compared with that of the oscillator bath, whose dissipa- 
tion rate is temperature- independent (for linear quantum 
Brownian motion). Indeed, the tanh (f2n/2fcBT) temper- 
ature dependence of the spin-bath dissipation rate (|26|) 
has been explicitly observed in systems such as glasses 
[32| where dissipation is mainly caused by interactions 
between phonon modes and TLS [33j . 

A similar argument also allows us to understand the 
absence of any temperature dependence of the normal- 
diffusion coefficient D [see Eq. (|27p] for the spin bath 
and thus of the rate of spatial decoherence (Fig. 
Eq. (|28b[) shows that for quantum Brownian motion 
this rate increases with temperature as coth(f2n/2fcBT) 
(and linearly with T in the high-temperature limit of 
the Caldeira-Leggett model [3]). This increase is due 
to the fact that, as the temperature is raised, increas- 
ingly excited energy levels will be occupied in each har- 
monic oscillator, and thus the characteristic wavelengths 
present in the bath will decrease. Shorter environmen- 
tal wavelengths mean that the bath will be able to 
better resolve the position of the central system, lead- 
ing to stronger decoherence (localization) of superposi- 
tions of well-separated positions. The quick saturation 
of the spin bath with increasing temperature implies 
that the characteristic wavelengths cannot become signif- 
icantly shorter, resulting for our model in a temperature- 
independent expression for D. 

The fact that the spin-bath decoherence rate D has a 
constant component implies that there is a heating source 
in the oscillator-spin model. This term in the master 
equation drives a diffusion process in the momentum and 
thus causes the average kinetic energy of the oscillator to 
increase linearly in time at a rate of at least Dq. This 
is likely to be a problematic source of heating in micron- 
scale ion traps. As this is independent of temperature, 
it will be apparent even if no systematic dissipation is 
observable (7 » 0). 



III. HARMONIC OSCILLATOR COUPLED TO 
A SINGLE TWO-LEVEL SYSTEM 
INTERACTING WITH A BOSONIC BATH 



F. Comparison between spin and oscillator baths 

In Fig. [5] we have plotted the temperature dependencies 
of the coefficients 7 and D for the spin bath [see Eqs. (j2"rj|) 
and (|27p ] and for an oscillator bath (quantum Brownian 
motion) [see Eqs. (|28a[) and (|28b[> ] . with both baths de- 
scribed by the spectral density (|2l>)) . Fig. [3k. shows that 
the spin-bath dissipation rate 7, Eq. (|2"rj)) , decreases with 



The discussion in the previous section indicates that 
the temperature dependence of the damping and dissi- 
pation rates arise from the fact that the bath seen by 
the central oscillator cannot absorb an arbitrarily large 
amount of energy. We are thus led to consider an ex- 
treme case in which the oscillator is coupled to only a 
single TLS which is itself weakly coupled to a bosonic 
bath. As mentioned in the Introduction, both in micron- 
scale, cryogenic ion traps and in GHz QEMS (the fre- 
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quency regime relevant to the observation of quantum 
effects) it is likely that only a few TLS take part in the 
dynamics, rather than a TLS environment with contin- 
uous spectrum as used in [2?], [28| . This lends particular 
experimental relevance to the single-TLS model consid- 
ered here. For example, our model could represent a 
single charge trap electrostatically coupled to a nanome- 
chanical resonator, with the charge trap itself coupled 
to Johnson-Nyquist electrical noise in the surrounding 
circuit. 

We begin with a single vibrational degree of freedom 
(the central harmonic oscillator S) coupled to a single 
spin- i -particle (TLS), 



H = Mloa^d + gX ® a z 



where 



-1/2 



X 



2MQ 0/ 

(2hMn y 1/2 P 



(a + at) = x, 
—i(a — a') = p, 



(29) 

(30) 
(31) 



are dimensionless position and momentum operators. In 
the case of a realization of the model as a charge trap, 
a z refers to two distinct charge configurations of a single 
microscopic trap in one of the electrodes near a nanome- 
chanical resonator or trapped ion. We could think of this 
as some kind of double well, in which case the eigenstates 
of a z refer to states localized on one side of the barrier or 
the other. The tunnel-split ground states under the bar- 
rier are eigenstates of a x and the tunnel splitting is A. If 
we define the eigenstates of a x by <r x |±) = ±|±), then we 
can define the "bit-flip" operator as a z = |+)(— 1 + |— )(+|- 
We expect that the TLS remains close to thermody- 
namic equilibrium with a heat bath at temperature T 
even in the presence of the coupling to the central oscil- 
lator. Thus we assume that its state at all times can be 
approximated by the thermal state 



where 



Ar=J>+|+)(+|+P-|->H J 



P± = c -A/fc B T^ 
P- 



(32) 



(33) 



The dynamical process that maintains the TLS in ther- 
modynamic equilibrium could be quite obscure. How- 
ever, a simple model can be given by weakly coupling 
the single TLS to a bosonic bath at temperature T. The 
coupling is capacitive, 



-^coupling = Oz ® ^2 9kqk{t), 
k 



(34) 



while the free Hamiltonian for the bath is a sum over 
harmonic oscillators, each with canonical coordinates 
QkiPk and frequency uif.. If the coupling is weak so that 
gu <C A, Wfc, the corresponding Markov master equation 



for the density operator p of the joint oscillator-TLS sys- 
tem is Hi 



= -ifio [a*a,p] - iA [a x ,p] - ig[Xa z ,p] 
+ 7(n + \)V[&-]p + -fnV[a + ]p, 



(35) 



where f2o is the vibrational frequency of the central os- 
cillator and 7 determines the heating rate. The super- 
operator T>[A] is defined by 

V[A]p = ApA^--{A^Ap + pA^A), (36) 

and <t+ = ar_ = |+)(— | are raising and lowering opera- 
tors in the eigenstates of a x , i.e., the energy eigenstates. 
Finally, the parameter n is defined by 



(37) 



It is easy to see that in the absence of the coupling to 
the oscillator, the steady state for the TLS implied by 
Eq. (j35|) is simply given by Eq. ([33]) . 

We can now calculate a heating rate for the oscillator. 
We do this by adiabatic elimination [35[ of the TLS, i.e., 
we assume that 7 3> flo, g, A so that the TLS remains in 
steady state slaved to the motion of the oscillator. We 
define the operators acting on the vibrational degree of 
freedom by p ++ = (+\p\+), p — = {-\p\~), p+- 

■ i) Note that the reduced density operator for the 
vibrational degree of freedom is just given by ps = p++ + 
p The master equation (|35p then implies that 



dp 



++ 



ig( x P-+ - P+- x ) - ifio [6*6, P++] 
- 7(n + l)p++ + jnp — , (38a) 

= -ig(Xp+- - P-+X) - in [ata, p—] 
+ 7(n + l)p++ - jnp — , (38b) 



dt 



dt 



-ig( x P — - P++ x ) - itio[ata, p^ 
-|(2n + l) / 5+_-2iAp + „. 



(38c) 



Assuming that the off-diagonal operators p^ reach a 

steady state, we find that 



P+- 



2ig 



7(2n+ 1 



-(Xp__-p ++ X). 



(39) 



Substituting this into the equation of motion for the di- 
agonal components leads to 

^- = -iCl [tfa,p s ] -T[X, [X,ps]], (40) 

where the last term implies diffusive heating and the rate 
is given by 



2g 2 



7(2n + l)' 



(41) 



7 



(a) 7 = 10 




(b) 7 = 100 




FIG. 3: Comparison of the adiabatic approximation (crosses) 
and the full master equation (circles) for the simple model 
discussed in Sec. IIIII We show the mean occupation number 
(N) of the central harmonic oscillator as a function of time at 
three different temperatures (a larger value of n, see Eq. (|37[1 , 
corresponds to higher temperature). We use A = 1, flo = 1, 
g = l, and (a) 7 = 10, (b) 7 = 100. 



a z , the vibrational degree of freedom suffers momentum 
kicks of equal magnitude but random sign. As the tem- 
perature goes to infinity, the TLS state is the identity 
operator and the vibrational degree of freedom suffers no 
kicks at all. For this reason the heating rate goes to zero 
at high temperature — a reappearance of the damping- 
rate feature of the oscillator-spin model discussed in 
Sec. Ill Fl Also note that at zero temperature there is 
a fixed momentum-diffusion rate that causes the oscilla- 
tor to heat, which again is the case for the full spin-bath 
model of Sec. |TTJ 



IV. SUMMARY AND CONCLUSIONS 

We have derived the master equations for a single har- 
monic oscillator coupled (i) to a bath of two-level systems 
(the oscillator-spin model) and (ii) to a single two-level 
system interacting with a bosonic bath. These mod- 
els and the derivation of the relevant master equations 
not only close an important gap in the set of canonical 
system-environment models for decoherence and dissipa- 
tion, but are also motivated by and relevant to current 
experiments such as quantum-electromechanical systems 
and micron-scale ion traps. 

For both models the key features that arise are: (i) The 
systematic dissipation rate for the oscillator decreases 
with increasing temperature; (ii) at zero temperature the 
systematic dissipation rate goes to zero; but (iii) there is 
a temperature-independent momentum-diffusion rate (a 
heating rate). Interestingly, this behavior is very differ- 
ent from the model for quantum Brownian motion, and 
we have explained how it can be understood as arising 
from a rapid saturation of the spin environment. 

An obvious direction for future investigations is the 
application of our model to concrete experiments and 
a comparison of the theoretical predictions with experi- 
mental data. In particular, we plan to carry out direct 
simulations of the model in an ion trap. 



A full numerical simulation (see Fig. [3|) confirms that this 
master equation is a good description of the dynamics in 
the limit of large 7. However, even when the adiabatic 
approximation is not the same as the full master equa- 
tion, the full model still shows diffusive heating. 

At zero temperature the TLS is in the ground state 
— ). As this is a superpositon of the two eigenstates of 
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